Specific Heat Discontinuity in Impure Two-Band Superconductors 



O 
o 

(N 



S3 
O 
O 

Oh' 
33 

03 



O 

> 
m 

On 
O 
(N 
O 
-i— > 



i 

C3 
O 

o 



X 



Todor M. Mishonov, 1,2, Q Evgeni S. Penev, 2 Joseph 0. Indekeu^'Q and Valery L. Pokrovsky 3 '0 

Laboratorium voor Vaste-Stoffysica en Magnetisme, Katholieke Universiteit Leuven, 
Celestijnenlaan 200 D, B-3001 Leuven, Belgium 
® Faculty of Physics, Sofia University "St. Kliment Ohridski", 5 J. Bourchier Blvd., 1164 Sofia, Bulgaria 
Department of Physics, Texas A&M University, College Station, Texas 77843-4242 
and Landau Institute for Theoretical Physics, Kosygin Street 2, Russia, Moscow 117940 

(Dated: February 1, 2008) 

The Ginzburg-Landau coefficients, and the jump of the specific heat are calculated for a disordered 
two-band superconductor. We start with the analysis of a more general case arbitrary anisotropy. 
While the specific heat discontinuity at the critical temperature T c decreases with increasing disor- 
der, its ratio to the normal state specific heat at T c increases and slowly converges to the isotropic 
value. For a strong disorder the deviation from the isotropic value is proportional to the elastic 
electron scattering time. In the case of a two-band superconductor we apply a simplified model 
of the interaction independent on momentum within a band. In the framework of this model all 
thermodynamic values can be found explicitly at any value of the scattering rate. This solution 
explains the sample dependence of the specific heat discontinuity in MgB2 and the influence of the 
disorder on the critical temperature. 

PACS numbers: 74.20.Fg, 74.72.-h 



I. INTRODUCTION 

The investigation of the specific heat C(T) is an im- 
portant tool for understanding the nature of the su- 
perconductivity and anisotropy of the superconducting 
gap A p (T) on the Fermi surface e p — Ep. Historically 
the relative specific heat jump AC / C n {T c ) was used to 
establish the BCS picture 1 for the conventional super- 
conductors having nearly isotropic gap. Subsequently 
the thermodynamics of clean anisotropic-gap supercon- 
ductors was analyzed in the weak coupling approxima- 
tion by Pokrovsky and Ryvkini^ They have found that 
anisotropy suppresses the value AC / C n (T c ) in compari- 
son to its isotropic value 1.43. This inequality is not satis- 
fied in classical low-temperature superconductors partly 
because they are not extremely clean, but also since the 
weak coupling approximation has a poor precision. Gei- 
likman and Kresini^ have proved that the first correc- 
tion due to interaction increases AC/C n (T c ) and thus 
disguises the anisotropy effect. The modern supercon- 
ductors display really high anisotropy. In particular, su- 
perconductivity is highly anisotropic in MgB 2 . This fact 
is the main motivation of this work. It is well known 
that only superconducting crystals of very high quality 
can reach the theoretical clean-limit asymptotics. As a 
rule, the reduced specific heat jump is sample-dependent 
and understanding of this disorder dependence is a chal- 
lenging problem. The latter is especially important for 
MgB2, a compound now being in the limelight of super- 
conductor materials science. 

The aim of the present paper is to derive the depen- 
dence of the critical temperature and the relative specific 
heat jump AC/C n (T c ) on the elastic scattering time of 
the charge carriers at the critical temperature r(T c ), for 



two-band superconductors having in mind application to 
MgB2. For this purpose we need corresponding formu- 
lae for a general dirty anisotropic superconductor. Such 
equations were derived inl We reproduce them here for 
readers convenience and because there occurred several 
misprints in the cited work which we correct here. In the 
r-approximation the electrical resistivity of the normal 
metal p c i is determined by this scattering time. Thus, 
our formula can be used for the investigation of cor- 
relations in the experimentally determined AC / C n (T c ) 
versus p e \(T c )/T c plot. The comparison of the theoreti- 
cal curve and the experimental data can reveal the gap 
anisotropy Ab yP and the scattering rate l/r(T c ). The gap 
can depend on both the quasimomentum p and the band 
index b. 



The applicability of the weak-coupling theory to MgpJ2 
is contentious. However, experimental results on the 
relative specific heat discontinuity indicate that the 
anisotropy effect is more profound than the effect of in- 
teraction. For this compound, the reduced specific heat 
AC/C n (T c ) is definitely smaller than the weak-coupling 
BCS value 1.43. Moreover the temperature dependence 
of the specific heat of the superconducting phase is de- 
scribed fairly well& by the two-band modeli and the rel- 
ative specific heat jump agrees^ with the Moskalenko's 
weak-coupling formula. 7 The comparison of the latter 
and the ab initio strong-coupling calculations 8 for MgB2 
shows that the decrease of AC/C n (T c ) due to different 
values of the superconducting gap for different bands is 
at least 2 times bigger than the increase of this reduced 
specific heat jump due to the strong coupling effects. We 
discuss this point in the concluding discussion. 



2 



II. CLEAN SUPERCONDUCTORS 

In this section we reproduce some results for an 
anisotropic clean superconductors obtained by different 
authors many years ago and derive equation for the spe- 
cific jump discontinuity in this case. Though neither of 
these results is new, they are necessary for understand- 
ing nest sections.lt was shown in Ref. Ijflthat, within the 
framework of the weak-coupling theory, the order param- 
eter possesses the property of separability: 



A P (T) = H(T) Xp . 



(1) 



According to Eq. QJ, the temperature dependence char- 
acterized by the factor S(T) is separated from the angular 
dependence described by the factor Xp- The Ginzburg- 
Landau (GL) expansion for the free energy density 910-11 
can be written in terms of the temperature-dependent 
factor 3(T) alone: 



f(B,T) = a 



T — T 1 
^|HP + i»|E|*. 



The specific heat jump per unit volume is related to the 
GL coefficients by the following relation: 



(C s ~C n )\ 



AC 



T c b 



(3) 



where C s is the specific heat per unit volume of the su- 
perconducting phase and C n is that of the normal phase. 

Our starting point are the expressions of Gor'kov and 
Melik-Barkhudaro\*i£ for the GL coefficients in the clean 
limit which can be written as 



/ 2\ - C(3, / 4\ 



where the Hurwitz and the Riemann zeta functions, 
Q{k,z) and C(^)> respectively, read 



C(M) = $> + *)-*, C(*) = C(M) = X/ 



-k 



(5) 



n=0 



n=l 



and obey the relation £ (k, |) = (2 fe — l)C(k). A simple 
variational derivation of Eq. (J1J is given in Ref. El The 
celebrity of £(3) in mathematics has been discussed in 
Ref. El The normalized moments of the gap anisotropy 
function are determined by averaging over the Fermi sur- 
face, having the general form in the D-dimensional case 



BZ 



e„=E F 



(6) 



space integral is taken over the whole Brillouin zone (BZ). 
The integration over the Fermi surface e p = E-p implicitly 
includes summation over fragments and sheets of differ- 
ent bands, if any. The normalizing factor 



(7) 



is the density of states (DOS) per unit volume for fixed 
spin, and enters the normal-phase specific heat 



C n (T) = -n 2 k 2 B v F T. 



(8) 



This equation together with the formulae for the GL co- 
efficients, Eq. (gj, lead to the following value for the re- 
duced jump of the specific heat 



AC 

(2) 7^ 



r i _(a^) 2 _(x 2 ) 

_L Pa (A*) 

7C(3)/3 a ' I 12 



12 



(X 4 ) 



< 1, 



I 7C(3) 



1.42613. 



(9) 

which is exactly the result obtained in Refs. 21:3; for a 
methodical derivation see Ref. El Using Eq. © for T 
slightly lower than T c , we get for the equilibrium order 
parameter 



ISI 2 



T -T c a 



|A r 



1 1—1 leqXp 



t c b ' 
2(27rfc B r c ) 2 T c -r (x 2 ) 2 



(4) which is the result by Gor'kov and Melik-Barkhudarovii£ 



III. DISORDERED ANISOTROPIC 
SUPERCONDUCTORS 



A. Transition line and order parameter 

In this subsection we analyze the transition tempera- 
ture T c as a function of the elastic scattering rate 1/r 
and the angular dependence of the order parameter \p- 
As it was explained before, the angular dependence is 
the same for any temperature at fixed r. The transition 
line has been studied in Ref. Although the equations 
obtained in the latter work were rather general, their 
treatment was focused on a specific situation — a mixture 
of s- and d-pairing characteristic for cuprate supercon- 
ductors. Therefore, it is useful to analyze the results for 
a less exotic case of anisotropic s-pairing. The general 
equation for the transition line found in Ref. reads 



where dS p is an infinitesimal surface element and v p = 
V p £p is the quasiparticle velocity. The quasi-momentum 
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where following notations are introduced: 



f(T,r) = - 



9(T,t) = - 

7T 



In- 



2wk H T 
I I ■'■ + \ 



"I 



(12) 
(13) 



x = (2nk B TT/h)^ 1 ; e is the cutoff energy; F(x) is the 
Euler digamma function; V n are eigenvalues of the linear 
operator V with kernel V(p, p') equal to the electron- 
electron effective interaction energy at the Fermi sur- 
face multiplied by the DOS z/p; v E'n(p) are the cor- 
responding eigenfunctions normalized according to the 
condition (I'I'nl 2 ) = 1. The transition temperature of 
the clean superconductor is determined by the equation 
f{T c0 , t — oo) = Vq -1 which gives 

k B T c0 = ^ exp (-n/Vo) , x = *o, (14) 

7T 



where 



7 = exp C 



4exp(F(l/2)) 



1.781, (15) 



C = 0.577 is the Euler constant, and Vq is the max- 
imum eigenvalue of the operator V. The angular de- 
pendence of the order parameter in this case is given 
by the corresponding eigenfunction 'i'o(p)- As long as 
x = (27rfc B T c0 T/S) -1 remains large the transition tem- 
perature T c differs from T c q by insignificant corrections 
of the order of Xq . We call the superconductor moder- 
ately dirty if the value x becomes small, but /(T c ,r) 
is still close to V . More precisely, it means that 
lner/ft - ttV^ 1 | < Vq 1 . Then the solution of Eq. JTTJ 
reads 



(16) 



T c (r) = T c0 - 



where k — (x 2 ) I (x) 2 > 1 is the anisotropy coefficient, 
which is equal to 1 for an isotropic superconductor, and 
To = 2hrf /(jrk B T c o). In this range of the scattering rate 
the angular dependence of the gap %p is identical to that 
for the clean superconductor: Xp = ^o(p)< Thus, the 
transition temperature decreases in a power-like way with 
the increase of the scattering rate 1 /r or the residual re- 
sistivity p ICS proportional to this rate. This is a pecu- 
liarity of the anisotropic superconductor. The exponent 
in Eq. 1)16(1 is zero for the isotropic superconductor, cf. 
Ref. |20j. Equation (|16|l was first derived by Hohenberg 21 
for weakly anisotropic superconductors. Its validity for 
arbitrary k in the range of moderate dirt was proved in 
Ref. 0. 

We call the dirt strong if the parameter In It /ttH be- 
comes less than V'o" 1 and has the order of magnitude of 
V , and if the difference V^ 1 — In er/irK is not small in 
comparison to Vq~ . Equation l|16|l remains qualitatively 
correct, but k becomes a slowly varying function of r. 



The exact formula for the transition temperature in this 
range is given by Eq. (37) of Ref. 0- 

In the extra-dirty limit -Iner/H becomes much 
smaller than V^ _1 , but still er/h^> 1. The last inequality 
ensures that the elastic scattering does not destroy the 
Fermi surface. In the extra-dirty limit the angular depen- 
dence of the gap reaches its limiting value Xp ^ ^(p)j 
where V(p) = (V(p, p')) ,. The equation for T c in the 
extra-dirty limit reads: 



k B T c (r) = — exp 

7T 



(V) 



(er/ny 



(17) 



where (V) = \V(p)) . It should be noted that k in the 
last equation differs from that for the clean superconduc- 
tor, namely: 



(EnVni^n)) 2 



(18) 



B. Specific heat discontinuity 

The theory of dirty anisotropic superconductors^ was 
based on the Green's functions method combined with 
the Abrikosov-Gor'kov averaging over the random impu- 
rity field^SiS A simplifying assumption was the isotropy 
of the scattering which is characterized by a constant rate 
1/t. In particular, the authors derived the GL equations 
and GL coefficients with an accuracy of a common scaling 
factor— (see also an earlier articled For the representation 
accepted in this article this factor is (x 2 )j as it follows 
from the comparison of Eq. 10} and Eqs. (48), (60), (78V- 
(82) in Ref.|j. Correcting a misprint in Eq. (59) in Ref. 4, 
further repeated in Eqs. (61), (82) therein, and slightly 
regrouping terms we find: 



(X 2 ) - Ux 2 ) - (X) 2 ) -cC: 



3 2,0 



(19) 



(7rfc B T c ) 



(x 4 )C3,o ~ (x 2 ) 2 a;cC4,o 
+4<X 3 )(X) XcUi 

+ 2 <X 2 )<X) 2 (^4,2+^3,2) 



where 



+ (x) 4 ^cC4,3 



K/T(T C ) Xq 



(20) 



2irk B T c T c jT C Q 



(21) 



is the dimensionless scattering rate extrapolated to the 
critical temperature. The resistivity of the normal metal 
Pol is determined by the Drude formula: 



Pel 1 



m 1 n tot e 2 T, 



(22) 
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where m is the effective mass tensor and ntot is the den- 
sity of normal charge carriers; for clean crystals the total 
volume density of all charge carriers, electrons and holes, 
eritot — e(n e — rih), can be determined by the Hall con- 
stant TZh — l/en to t in strong magnetic fields. For clean 
superconductors, disregarding some subtleties, the same 
ratio n/m enters the London penetration depth A c i can at 
T = 



1 



lean 



(0) 



mc 2 £o 



(23) 



where in Gaussian units Eq = 1/4tt; A c i ean (0) ~ 0.1-1 /im. 
Multiplying these equations we obtain a useful estimate 



hc 2 Sop c \(T c ) 
2^A c 2 lcan (0)fc B TV 



(24) 



The notation Qk,l in Eq. fflfy stands for the generalized 
zeta- functions defined in Ref. and taken at the value of 
the argument z c = x c + 1/2, i.e^ 



Ck,l = Ckd(z c )- 



(25) 



For the readers' convenience we recall the definition of 
these functions: 



Ck,j(*) = X> + 1/2) 



-/ 



(26) 



n=0 



They represent a natural generalization of the Hurwitz 
zeta functions: 



Cfc, («) = c(M> a,o(i) = c(*,i) = c(fe). 



(27) 



Below we provide the asymptotics of Ck,i{z) for z — > oo 
necessary for the further calculations: 

Gw(*) - i)C(0*~ fc - fe[(2 ,_1 - i)C(i - 1) 

-(2'-lX(0/2]« _fc_1 , for ?>2, 

Cm( z ) ~-2~ fc Inz, for fc > 1, 

Cfc,o(*)~(fc-l)~ 1 *~ fc+1 . for fc>l. 

(28) 

Let us note that, for integer arguments k, the Hurwitz 
zeta functions are associated with the Euler polygamma 
function F^: 



C(fc + l,z) 



fc! 



F (fc) (z), fc = 1,2,3,. 



(29) 



With these notations the reduced discontinuity of the 
specific heat reads: 



AC 



k B T c b 



12 1 

703)^7' 



(30) 



where 



j r =7C(3) [<X 2 ) ~ ((X 2 ) - (X) 2 ) ^2,0 



(X 4 ) C3,0 - (X 2 ) 2 XcCifl 

+ 4<x 3 )(x) ^cC3,i 

+ 2( X 2 >(X) 2 (^C4,2+^C3,2) 



+ (X> 4 2cC4,3 



(31) 



This general equation will be applied in the following 
subsections to some important special cases. 

It should be stressed that for isotropic superconduc- 
tors, (x™) = 1, the specific heat jump is impurity inde- 
pendent: [3 T = 1. The proof is straightforward taking 
into account the identity 

C3,0 ~ Xc(4,0 + 4X C C 3 ,1 + 2 (x 3 c (^ 2 + X 2 c(3,2) + 2^4,3 

= C(3,1/2) = 7C(3). (32) 

Likewise, using Eq. I|28|l one can prove that the asymp- 
totic form of AC / C n {T c ) for an extremely disordered su- 
perconductor with an arbitrary anisotropy is given by, to 
leading order in x~ , 



AC 



12 



C n (T c ) 7C(3) 



(*-l) ( 



5tH 



x c \14C(3) 



(33) 



C. Two-band superconductors 

1. Critical curve and order parameter 

Keeping in mind the application to MgB 2 (for a review 
see Ref. Ilq) . we apply the general results of the previous 
sections to a simplified model of a two-band supercon- 
ductor. In this model we assume that the Fermi surface 
consists of two disconnected sheets having different DOS. 
The interaction amplitude V(p, p') is assumed to be a 
constant within each band. Thus, it can be described by 
a 2 x 2 matrix: 



V = 



Wi u 
u w 2 



(34) 



where W\ , Wi are the interaction energies between any 
two points within the first and the second sheet of the 
Fermi surface, respectively; U is the interaction between 
any two points of different bands. 

Let us first work out the transition temperature 
T c q and the order parameter x f° r a clean two-band 
superconductor^ For our simplified model the order pa- 
rameter x(p) is a constant within each band, i.e. it can 
be represented by a 2-component vector: 



X 



5 



The eigenvectors of the operator V obey the following 
linear equations: 



Cif/*i + C 2 W 2 ^ 2 = A*2, 



(35) 



where the coefficients c X2 — V\^j{y\ + v 2 ) are the statis- 
tical weights of the two bands, which reflect the integral 
character of the operator V. The two independent eigen- 
values of Eqs. i|3"5|) read 



Vb,i = A± = rj±e, 



(36) 



where 



V = \{cxW x + c 2 W 2 ), e = Ve + cic 2 U 2 , 
£ = \{ciWx -c 2 W 2 ). 
The corresponding eigenvectors are: 



v sign(f/) v / 5 i 7 (l-e/e)/ ' 
-sign^^a-e/e)' 



(37) 



(38) 



2c 2 



This is apparently a kind of the Bogolyubov transforma- 
tion. Both vectors are normalized: 

(|*±| 2 ) = ci |*±i| 2 + c 2 |* ±2 | 2 = f . (39) 

The average values of the anisotropic eigenfunctions read: 



(40) 



= - sig n(^/| /f (l + f)- (41. 



ft is useful to write simple expressions for the squared 
averages: 

<*±} 2 ^±^[(ci-c 2 )£ + 2 ClC2 [/], (42) 
<* + ) 2 + <*_) 2 = l, { X f = (y + f. 
For the gap ratio 5 Eq. gives 

^^ = ii± = ^=sign(C/)^/i±|. (43) 

Whence the moments of the anisotropy function read 
(A") Cl S n + c 2 



(x n ) 



(A 2 )™ /2 M 2 + c 2 )™/2 



(44) 



Using the general results formulated earlier, we find the 
transition temperature of the clean two-band supercon- 
ductor: 



t c o = e 

7T 



(45) 



Equation 1)11(1 for the critical curve within this model can 
be simplified to the form: 

^9 [r) + (ci - c 2 )£ + 2cxc 2 f7 - c x c 2 fd\ 

= (1-A+/)(1-A_/). (46) 

Here we have denoted d — detV — W X W 2 — C/ 2 , 
and abbreviated the functions /(T c , r) and g{T c , r) as 
/ and g, respectively, cf. Eqs. Q12JI. ((13J1 . It is con- 
venient to introduce the dimcnsionless scattering rate 
xq = (2Trk B T C QT/H)~ 1 and the dimcnsionless transition 
temperature 9 — T c (t)/T c0 ; x c — xq/6. In terms of these 
variables the functions / and g read: 

f = ^---\n9-g, 



irg = F 



I 



(47) 



The equation for T c finally reads 



(x) 2 A+ + (l-(x> 2 ) A_-/A + A_ 

= (1-/A+)(1-/A_), (48) 



where 



2 (ci<5 + c 2 ) 
C\0 A + c 2 



C2 



1-Cl. 



(49) 



To give a realization of the possible dependence 9(x ), 
we have made a numerical calculation setting^ S = 2.63, 
A + = 1.02, A_ = 0.45, and c x = 0.42 2; oth er authors give 
slightly different values, cf. Refs. l5l l6ll7lll8lll9L The re- 
sults for 0(xq) and AC / C n (T c ) vs x c are shown in Fig.^ 
In the asymptotic regions of moderate and extreme dirt 
Eq. ijTT)|) is valid with 



I" 



(ci - c 2 )£ 



cic 2 C/ 



Cl 5 2 



C2 



(ci<5 + c 2 )^ 



for the clean and moderate dirt cases, and 

c\Wl + c^VF 2 + c x c 2 (2r) + £/)£/ 
(c 2 VFi+c 2 IF 2 + 2c lC2 [/) 2 

for the extreme dirt case. The order parameter in the 
moderate dirt range is proportional to , J , + . In the range 
of strong disorder it reads: 



c\W\ + c 2 U - it x c\c%d ln(er//i) 
c\U + c 2 W 2 — 7r _1 cic 2 (i \n(eT/h) 



(50) 



6 
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FIG. 1: (a) Critical curve T c (t) for two-band model: the re- 
duced critical temperature 6 — T c (t)/T c q versus dimension- 
less scattering rate xo = h/2nk B T c oT oc p c \(T c ). The set of 
parameters corresponds to MgB2; for details see the text, (b) 
Reduced specific heat jump AC/C n (T c ) as a function of the 
dimensionless scattering rate x c = h/2nk B T c r oc p c \(T c )/T c 
for the same set of parameters. The dashed line indicates the 
asymptotic (BCS) value for x c S> 1. The curve will be shifted 
up by about 20% due to the strong-coupling effects. 



In the limit of extreme disorder \n(eT/h) <C W 1 \ it tends 
to the limiting value: 



Xoxtr = V 



c{W\ + c 2 U 
c 1 U + c 2 W 2 



(51) 



It is worthwhile to note that at W\ = W 2 = U the 
anisotropy parameter k is equal to 1 for any value of the 
scattering rate independently on the values C\ , c 2 , and all 
thermodynamical values do not depend on r similarly to 
the completely isotropic case. 



2. Specific heat discontinuity 

We note that the normalization factors of the gap- 
anisotropy function and the superconducting gap cancel 
each other in the formulae for the experimentally measur- 
able jump of the specific heat. Therefore, we can use the 
normalization (x 2 ) — 1 without loss of generality. The 
calculations made above can be used directly to find the 
averages necessary for the calculation of the specific heat. 
The influence of the disorder on the relative specific heat 
discontinuity for a two-band superconductor is shown in 
Fig. n (b) ■ In the extreme dirty limit we find that the 
relative specific heat discontinuity tends to its isotropic 
value in agreement with the fact that the density of states 
becomes isotropic in this limit^ 



D. Dirty isotropic alloys 

It is straightforward to verify that in the isotropic case 
(x = const on the Fermi surface) the coefficients ao and 
b do not depend on r. Thus, neither the energy gap, nor 
the specific heat are influenced by impurities, in accor- 
dance with the Anderson theorem^ Analyzing Eq. 
we conclude that the transition temperature also does not 
depend on the scattering rate. Indeed, in the isotropic 
case the eigenfunctions of the operator V are spherical 
harmonics = Yi m {9, ip). Among the latter only one, 
^007 has a nonzero average. Thus, Eq. takes the 
simple form: 



g(T c ,T) + f(T c ,r) = V Q 



(52) 



According to the definitions (|12ll and l|13l) . the sum 

/(T,t) + g{T,T) = ln^^f^ - F(l/2) does not depend 

on the scattering rate. Hence, T c also does not depend 
on the scattering rate. 



E. Separable approximation 

The separable approximation 

^(P,P') = £V«*n(P)*n(p') « VoXpXp', (53) 

n 

where ^E'o(p) = Xp is very often used for modelling of 
the gap anisotropy in superconductors A p « S(T)x P - 
As we demonstrated earlier, this approximation is valid 
in the range of clean and moderately dirty superconduc- 
tors. Applying this approximation to the equation for T c 
in two-band model Eq. (|48|l we obtain Moskalenko and 
Palistrant, 31 Abrikosov 32 and Kogan's equation^ 3 . 



In 



T, 



r0 



T e (r) 



= 1- 



(A P ) 2 
<A P ) 



(54) 
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where^i 1/r = §(1/ti2 + I/T21), and 1/ti 2 and l/r 2 i are 
rates of interband scattering. The results of the numeri- 
cal solution of this equation are depicted on Fig. [3 For 
superconductors with zero averaged gap (A p ) = which 
are p- and d-type superconductors, for example, this 
equation formally coincides with the Abrikosov-Gor'kov 
result 2 ^*^ for superconductors with magnetic impuri- 
ties; superconductivity disappears at the critical value 
xq = 1/47 = 0.1404. For weak disorder we have 



Tco — T c 



(X 2 )-(X) 2 TTfi 

(X 2 ) 4fc B r 



«T C , 



(55) 



In such a way one of the most important properties of 
multigap and anisotropic superconductors is that the 
nonmagnetic impurities are pair breaking, similar to 
magnetic impurities in conventional superconductors^^ 
A similar influence of structural defects was discussed 
by Abrikoso\sS for triplet superfluids. The reduction 
of the critical temperature by disorder has been ob- 
served for layered cuprateS(2Si2L2Si2S for impurity scat- 
tering in triplet superconductor— UPt3 and recently for 
MgB2. 16 ' Only dimensionless ratios of the gap function 

moments, like (x) 2 / (x 2 ) m Eq. J52J|, or (x 2 ) 2 / (x 4 ) 
in Eq. are relevant for the thermodynamics of su- 
perconductors. This explains why strongly anisotropic- 
gap layered cuprates were seemingly successfully ana- 
lyzed as two-band superconductors^ (this reference is 
a comprehensive review of the properties of multigap su- 
perconductors), and vice versa why the first prominent 
two-gap superconductor MgB2 could be analyzed as if it 
were a single-band anisotropic-gap superconductor^ cf. 
Ref . M 



IV. DISCUSSION AND SUMMARY 

In order to investigate AC/C n (T c ) versus 1/tT c cor- 
relations it is necessary to have a good method for the 
determination of r. This could be far-infrared measure- 
ment of the high-frequency conductivity, or merely the 
static resistivity. For polycrystalline materials one has 
to use compounds like MgB 2 with good connections be- 
tween grains. The variation of the residual resistivity can 
be produced by radiation. 

For MgB2 the gap ratio S = A1/A2 can be determined 
by spectroscopic measurements, and the ratio of the DOS 
V1/V2 can be calculated from first principles. Under these 
conditions, for high-quality clean samples we can evaluate 
the up-shift of the AC / C n (T c ) curve, and the influence of 
strong-coupling effects. Another complication is related 
to the variation of the electron wave functions in the two 
bands. This can lead to different scattering rates in the 
two bands. We consider that the scattering rate in the a- 
band with bigger gap is more important for the disorder 
reduction of the specific heat jump. If the scattering time 
cannot be determined by spectroscopic measurements we 
can use the value of the resistivity at the critical temper- 
ature p \(T c )/T c . Thus, the dependence of T c and the 
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FIG. 2: Critical curves 9 = T c /T c0 vs x = h/2nk B T c0 T for 
different gap anisotropies {xf / (x 2 ) = 0, 0.2, 0.4, 0.6, 0.8, 
and 1. For an isotropic superconductor, (x) = 1, the critical 
temperature is disorder independent while for p- or d-type su- 
perconductors (as some CuCb superconductors are) (x) = 
and the transition line is exactly the same as the Abrikosov- 
Gor'kov—*— curve for magnetic impurities in isotropic super- 
conductors. 



relative specific heat jump on resistivity at T — T c will 
be given by our formula and Fig-H Fig.Q](b), with the 
scale of the abscissa being a fitting parameter. We expect 
that the derived weak-coupling formula can be as useful 
for the analysis of AC versus p c \ correlations as the weak- 
coupling theory was successful in describing the temper- 
ature dependence of the specific heat for the clean MgB 2 
samples. 6 Thus, we conclude that the weak-coupling the- 
ory of the impurity reduction of the specific heat jump 
can reveal the main trend and qualitative properties of 
the effect. 

Irradiated superconductors are also a good example for 
the application of the present theory. Furthermore, let us 
note that conventional dirty superconducting alloys, for 
which a big enough series of samples with continuously 
changing resistivity can be prepared, are the best tools 
to investigate the influence of disorder on the thermody- 
namics of superconductors. 

Finally, let us summarize our results, (z) In anisotropic 
superconductors the transition temperature is suppressed 
by disorder like T c ~ r K_1 , where k = (x 2 ) I (x) 2 is an 
anisotropy parameter which is a slowly varying function 
of t. (ii) The order parameter retains its angular depen- 
dence until er 3> 1, whereas the DOS becomes isotropic 
for lner <c . The anisotropy of the order parame- 
ter can be probed in tunneling experiments. (Hi) The 
specific heat is suppressed by disorder just like T c , with 
an accuracy of a slowly varying factor, (iv) The rela- 
tive jump of the specific heat AC / C n (T c ) is smaller than 
its isotropic value in the clean limit. It is enhanced by 
disorder tending to its isotropic limit in the extreme dis- 
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ordered case, (v) In isotropic superconductors, and the 
two-band model with Ai = A 2 and arbitrary Cx/c2, T c 
and AC do not depend on the scattering rate 1/r; in 
particular this is the case for W\ — W% = U. 
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